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Abstract. The type II supernova is considered as a candidate site for the production 
, of heavy elements. Since the supernova produces an intense neutrino flux, neutrino 

OO ■ scattering processes will impact element formation. We examine active-sterile neutrino 

conversion in this environment and find that it may help to produce the requisite 
neutron-to-seed ratio for synthesis of the r-process elements. 



The r-process of nucleosynthesis accounts for the most neutron rich of the heavy 
elements. The most likely environment for this type of synthesis is the late time 
(t > 10 s post-core bounce) supernova environment. Many studies have explored 
this 'neutrino driven wind' as a candidate environment and found it to be poten- 
tially viable [1,2]. However, to date, no model correctly reproduces the observed 
abundance pattern. 

In the neutrino driven wind, material in the form of free nucleons is 'lifted' off of 
the surface of the neutron star by energy deposited by neutrino interactions. Ana- 
lytic and semianalytic parameterizations of the thermodynamic and hydrodynamic 
conditions in the wind can be obtained [3,4]. Models of this type may be used to 
explore the range of conditions within the context of the wind which will produce 
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the solar system distribution of r-process elements. The key determinant of whether 
a given scenario will produce the r-process is the neutron to seed nucleus ratio at 
the onset of the neutron capture phase. This ratio must be quite high (R > 100) 
in order to produce the very neutron-rich r-process elements. The factors which 
determine the neutron-to-seed ratio are the entropy of the material, the hydrody- 
namic outflow timescale and the electron fraction, Y e = 1/(1 + n/p) where n/p is 
the neutron-to-proton ratio. A study of many possible model parameters shows 
that one must decrease the electron fraction, and/or increase the entropy and/or 
decrease the hydrodynamic outflow timescale, relative to the conditions found in 
typical wind models, in order to produce the neutron-to-seed ratio necessary for 
the r-process [5,6]. 

Including the effects of neutrino interactions in general tends to make the requi- 
site conditions for r-process element production more extreme [7,8]. In particular, 
neutrino capture on free nucleons during alpha particle formation increases the elec- 



tron fraction [8]. This is the "alpha effect". Other neutrino process are discussed 
in [2,9]. 

There are three possible solutions to this problem. The first is that the supernova 
is the site of r-process synthesis, but it does not occur in the neutrino driven wind 
as it is currently modeled. The second is that the r-process elements are made 
at some other site such as neutron star-neutron star mergers. However, timescale 
arguments combined with isotopic abundance measurements and observations of 
old halo star metallicity show that this site is unlikely to account for the entire 
r-process distribution [10,11]. 

The third solution is the one that is investigated here: active-sterile [y e <-> u s , 
z7 e <-> z7 s ) neutrino transformation. The v s in our study is defined as a particle which 
mixes with the v e (and possibly also with and/or u T ) but does not contribute 
to the width of the Z boson. 

If we neglect v-v forward scattering contributions to the weak potentials, then 
the equation which governs the evolution of the neutrinos as they pass though the 
material in the wind can be written as: 
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The upper sign is relevant for neutrino transformations; the lower one is for antineu- 
trinos. In these equations 5m 2 = m\ — m\ is the vacuum mass-squared splitting, 
9 V is the vacuum mixing angle, Gf is the Fermi constant, and N~(r), N+(r), and 
N n (r) are the total proper number densities of electrons, positrons, and neutrons 
respectively in the medium. Resonances can occur when the on-diagonal terms in 
the wave equation are zero. The quantity in the brackets in Eq. 2 is proportional 
to Ye — 1/3. Since the electron fraction can take on values between zero and one, 
the bracketed quantity in Eq. 2 can be either positive or negative. Therefore, for a 
given choice of Sm 2 , resonances may occur for neutrinos or antineutrinos depending 
on the value of the electron fraction. 

In order to determine the survival probabilities of the neutrinos and antineutri- 
nos, we must know the electron fraction. In the neutrino driven wind, neutrino 
and antineutrino capture are the most important reactions in determining the elec- 
tron fraction. However, near the surface of the protoneutron star there is also a 
contribution from electrons and positrons: 



z/e + n^p + e ; v e + p # n + e 
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Therefore, the problem involves a feedback mechanism. The neutrino capture rates 
determine the electron fraction and the electron fraction determines the potential 



which is used in the neutrino transformation equations. These equations then 
determine survival probabilities of the neutrinos, and therefore their capture rates. 

We note that we can not assume weak equilibrium, (very fast capture rates 
compared with outflow rate) or weak freeze-out (very slow rates compared with the 
outflow rate). The rates must be tracked numerically. A previous study took the 
limit of weak equilibrium and found different behavior [12] than presented here. 

We perform calculations by tracking a mass element in the neutrino driven wind. 
We use the results of analytic models [4], where r oc exp(— t/r) and p oc r -3 
where r is the outflow timescale. Close to the surface, before the wind begins to 
operate we use the density profile of Wilson and Mayle [13]. At each time step 
we calculate all thermodynamic quantities, calculate the weak rates and evolve the 
neutrino transformation equations forward. We assume, since the outflow timescale 
is short, t ~ r < 0.5 s, that each mass element will have the same evolution as the 
previous one. Since we use a nuclear statistical equilibrium calculation, we cut 
our calculations off at the time when heavy nuclei begin to form. More detail is 
contained in [14]. 

A calculation for one mass element is shown in Fig. 1. The upper curve shows 
the evolution for the case of no neutrino oscillations. The wind parameters were 
r = 0.3 s, s = 100. The solid line shows the actual electron fraction, while the 



FIGURE 1. The electron fraction is plotted against distance from the center of the neutron star. 
The upper line shows the evolution with no transformation. The lower line shows the evolution 
of active-sterile mixing parameters of sin 2 9 = 0.01 and Sm 2 = 20 eV 2 . The dotted line shows the 
value that the electron fraction would take on if weak equilibrium obtained. 
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dotted line shows what the electron fraction would be if weak equilibrium obtained. 
The initial rise is due to Pauli unblocking of the electrons at the surface of the 
proto-neutron star. The lower curve shows the evolution of the electron fraction 
for mixing parameters of 5m 2 = 20 eV 2 , sin 2 2# v = 0.01. In the latter case there 
is a rapid drop in the electron fraction when the neutrinos begin to transform. 
In fact, there are three neutrino transformations. Initially, the electron neutrinos 
transform to steriles. Later antineutrinos transform, and finally, when the density 
falls far enough the antineutrinos transform back. The first transformation of the 
antineutrinos is seen in the small bump in the equilibrium electron fraction. A small 
"alpha effect" can be see as the slight rise in both solid lines at large distance. 

The active-sterile transformation scenario successfully reproduces a low electron 
fraction which is beneficial to the r-process. We now consider a range of 5m 2 , sin 2 9 
parameters. A contour plot of the electron fraction is shown in Fig. 2. Inside the 
dashed contour shows the region where conditions are neutron-rich enough to be 
favorable for r-process nucleosynthesis. In the bottom left corner of the plot, the 
solution is approaching the case without neutrino transformation. 

Although not shown here, we have studied a range of timescales for the neu- 
trino driven wind models, and seen that the qualitative features of this effect are 
reproduced [14]. 

We have used several approximations in our calculations, which we are contin- 
uing to study. These include the importance of the neutrino-neutrino scattering 




FIGURE 2. Contour plot of electron fraction as measured at the point where heavy nuclei begin 
to form. Neutrino driven wind parameters employed here are s/k = 100, r = 0.3s. . 



background in the oscillation equations, the importance of nonradial paths and 
feedback in the dense region near the proto-neutron star. The region where the 
latter two are important is the shaded region in Fig. 2. These problems are being 
studied in [15]. 

Conclusions: Meteoritic and observational evidence points to supernovae as the 
source of the r-process elements, although a self-consistent model of the neutrino 
driven wind which will produce the r-process elements is still elusive. In the next 
few years significant advances are expected in supernova modeling. If and when 
potential hydrodynamic solutions are exhausted and the caveats above have been 
explored, then the r-process may provide a signature for new neutrino physics. 
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